clc;clear;close all;
group = 5;
RMSE_PHI = zeros(group,4);
for g = 1:group
    Actan_RMSE_PHI = load(['.\Data\Data (',int2str(g),')\Actan\',int2str(g),'-RMSE_unwarpphase.mat']);
    RMSE_PHI(g,1) = Actan_RMSE_PHI.RMSE_unwarpphase;

    AIA_RMSE_PHI= load(['.\Data\Data (',int2str(g),')\AIA\',int2str(g),'-RMSE_unwarpphase(AIA).mat']);
    RMSE_PHI(g,2) = AIA_RMSE_PHI.RMSE_unwarpphase;

    TLS_AIA_RMSE_PHI = load(['.\Data\Data (',int2str(g),')\TLS_AIA\',int2str(g),'-RMSE_unwarpphase(TLS_AIA).mat']);
    RMSE_PHI(g,3) = TLS_AIA_RMSE_PHI.RMSE_unwarpphase;
%     
    TLS_LRD_RMSE_PHI = load(['.\Data\Data (',int2str(g),')\TLS_LRD\',int2str(g),'-RMSE_unwarpphase(TLS_LRD).mat']);
    RMSE_PHI(g,4) = TLS_LRD_RMSE_PHI.RMSE_unwarpphase;
    
end

x = linspace(1,group,group);
figure;hold on;grid on;
plot(x,RMSE_PHI(:,2),x,RMSE_PHI(:,3),x,RMSE_PHI(:,4));
legend('AIA','TLS','LRD');
hold off;


